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HIGHLIGHTS 


•  Physical  and  electrical  behavior  of  alkaline  electrolyser  is  simulated. 

•  An  electrical  analogy  is  used  to  emulate  the  behavior  of  the  alkaline  electrolyser. 

•  The  impact  of  the  operating  temperature  and  pressure  has  been  evaluated. 

•  The  effect  of  gas  bubbles  adhering  on  the  electrodes  has  been  analyzed. 
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This  paper  presents  an  alkaline  electrolyser  simulation  tool  (AEST).  This  new  approach  is  based  on  a 
physics  model  of  an  electrolytic  cell  and  an  electrical  analogy  of  the  physical  effects  taking  place  in  the 
electrolyser.  The  physics  model  has  been  developed  taking  into  account  the  electrodes’  electrochemical 
properties,  the  thermodynamics  of  the  reaction,  and  the  electrical  properties  of  the  electrode  plates, 
membrane  and  electrolyte.  The  electrical  analogy  is  used  to  integrate  the  physics  model  into  the  elec¬ 
trical  simulation  systems.  The  AEST  is  able  to  describe  the  electrolyser’s  physical  and  electrical  perfor¬ 
mance  and  estimates  the  operating  current  as  function  of  the  voltage,  temperature,  pressure  and 
electrolyte  concentration.  The  technical  information  on  the  Hydrogen  Research  Institute  (HRI)  electro¬ 
lyser  has  been  used  to  set  up  the  AEST.  The  start-up  phase  of  the  electrolyser  has  been  analyzed  for 
various  operating  temperatures,  pressures  and  voltages.  The  steady-state  response  of  the  electrolyser  (at 
constant  temperature,  pressure  and  concentration  of  electrolyte)  has  been  simulated  as  an  element  of  a 
possible  hybrid  system.  Simulations  have  been  carried  out  using  MATLAB/Simulink/SimPowerSystems®. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Nowadays,  renewable  energy  systems  combined  with  hydrogen 
production  are  being  considered  as  a  promising  source  of  clean 
energy,  ahead  of  fossil-based  sources.  In  hydrogen  production 
based  on  water  electrolysis,  there  are  three  principal  types  of 
electrolysers:  alkaline  electrolyte,  proton  exchange  membranes 
and  high  temperature  steam.  Alkaline  electrolysers  are  the  most 
mature,  best-known  and  most  used  technology  for  hydrogen 
production. 

Alkaline  electrolyser  models  can  also  be  divided  into  two 
groups.  The  first  group  is  known  as  the  characteristic  models  and 
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the  second  one  as  the  electrical  models.  The  majority  of  the  char¬ 
acteristic  models  are  based  on  the  semi-experimental  coefficients 
found  by  fitting  the  current-voltage  (i— V)  curves  obtained  from  the 
experimental  data.  Some  interesting  work  in  this  field  has  been 
reported  by  Vanhanen  et  al.  and  Ulleberg  [1,2].  These  models  have 
been  used  to  obtain  the  electrolyser’s  i—V  curve  and  the  rate  of 
hydrogen  production  [3,4].  In  this  kind  of  model  the  operating 
voltage  is  the  result  of  the  current  applied  as  input.  Physics 
modeling  has  been  used  recently  [5-7]  for  alkaline  electrolysers. 
This  physics  approach  is  focused  on  modeling  separately  the  effects 
taking  place  in  an  electrolysis  cell. 

Electrical  models  of  electrolysers  are  used  to  integrate  these  into 
electrical  simulation  systems.  These  electrical  representations  have 
been  proposed  in  different  research  publications  [8—11].  However 
most  of  those  propose  a  simple  resistance  plus  voltage  source  as  the 
electrical  equivalent  system.  Other  authors  use  transfer  functions  to 
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simulate  the  electrolysers  in  power  systems  and  to  evaluate 
hydrogen  production  strategies,  control  systems  [12,13]  and  grid 
frequency  regulation.  This  modeling  approach  is  only  valid  for 
electrolysers  working  at  constant  power  in  a  hydrogen  production 
plant. 

This  paper  proposes  an  Alkaline  Electrolyser  Simulation  Tool 
(AEST)  based  on  a  physics  model  with  an  electrical  analogy  in  order 
to  emulate  the  physical  and  electrical  behavior  of  the  alkaline 
electrolyser. 

This  paper  specifically  describes  in  detail  the  proposed  physics 
model  for  one  alkaline  cell.  The  validation  of  the  physics  model 
using  the  Hydrogen  Research  Institute  (HRI)  electrolyser  [3]  and 
Phoebus  electrolyser  [2]  is  also  presented.  The  electrolyser’s  elec¬ 
trical  analogy  based  on  the  physics  model  is  explained.  Finally,  the 
simulation  results  using  the  proposed  simulation  tool  for  electro¬ 
lyser  start-up  phase  and  steady-state  operation  are  also  presented. 


the  exchange  current  densities  for  anode  and  cathode.  In  the 
literature,  they  are  generally  expressed  in  mA  crrT2. 


b  = 


2.303 RT 
nFa 


(4) 


where  a  is  the  transfer  coefficient  which  can  be  defined  for  the 
anode  and  the  cathode  respectively  as  aa  and  ac. 


(5) 


where  Se ff  is  the  effective  electrode  surface  in  cm2,  and  is  defined  by 
(6)  as  a  function  of  the  fractional  electrode  coverage  (6)  of  the  gas 
bubbles  adhering  on  the  electrodes  [16].  In  this  paper,  the  same 
value  of  6  is  applied  for  both  electrodes.  The  value  of  6  can  be 
calculated  using  the  empirical  correlation  proposed  in  Ref.  [7]. 


2.  Physics  model 

2.1.  Electrolytic  cell  model 

The  proposed  physics  model  considers  separately  every  physical 
phenomenon  taking  place  in  the  electrolysis  process  in  an  alkaline 
electrolyser.  In  this  approach,  the  electrical  current  (z),  the  elec¬ 
trolyser’s  temperature  (T)  and  pressure  (P),  and  the  electrolyte’s 
molar  concentration  (m)  have  been  set  as  the  model’s  main  pa¬ 
rameters.  The  i—V  relation  is  then  defined  by  (1). 

^cell  =  Eth  +  ^act -c-d  +  ^act -a-d  +  (^c  +  ^a  +  Kele  +  ^mem  W 

(1) 

where  VCeii  is  the  alkaline  cell  voltage,  Eth  is  the  required  voltage  to 
start  the  electrolysis  reaction,  Vact-c-0  and  ^act-a-0  are  the  electrodes’ 
(cathode  and  anode)  activation  overvoltages,  Rc  and  Ra  are  the 
electrode’s  plate  resistances,  Re \e.£  is  the  electrolyte’s  resistance, 
Rmem  is  the  membrane’s  resistance,  and  i  is  the  operating  current. 

2.1.1.  Water  electrolysis ’  theoretical  voltage 

The  voltage  Eth(r,p,m)  has  been  defined  using  Nerst’s  equation 
with  a  T,  zrz  and  P  dependence  (2).  The  first  term  of  the  equation  is  a 
thermodynamic  approximation  of  the  Gibbs  energy  (Erev(T,i))  -  in 
this  model  the  correlation  proposed  in  Refs.  [7,14]  has  been  used. 
The  second  part  of  the  equation  includes  the  pressure  and  the 
electrolyte  concentration  impact  on  the  voltage  [7]. 

Eth (T,P,m)  =  £rev(T,l)  ~  [(^  “  Ph2o)2/Ph2o/Ph2o]  (2) 

where  R  =  8.315  J  K_1  mol-1  is  the  universal  gas  constant,  rz  =  2  is 
the  number  of  electrons  transferred  in  the  electrolysis  reaction,  and 
F  =  96,485  C  mor1  is  the  Faraday  constant.  The  partial  pressure  of 
pure  water  vapor  is  given  by  p*H20  as  well  as  the  electrolyte  partial 
pressure  (electrolyte  having  a  molarity  m )  by  pn20  [7-14]. 

2.1.2.  Electrochemical  electrodes’  activation  overvoltage 

Tafel’s  approximation  has  been  used  to  estimate  the  overvoltage 
due  to  the  electrochemical  activation  of  the  electrode  [7],  which  is 
defined  as  follows: 

Vact-a  -  ^al°g(jg — ^act_c  =  Mog(j^ — ^  (3) 

The  overvoltage  on  the  cathode  is  defined  as  Vact-c  and  the  one 
on  the  anode  as  Vact.a.  ha  and  bc  are  the  Tafel  slope  for  anode  and 
cathode  respectively,  they  can  be  computed  by  using  (4).  J  is  the 
current  density  of  the  electrodes  15]  given  by  (5).  Jo-a  and  Jo-c  are 


Self  =  5(1  -  6)  (6) 

where  S  is  the  nominal  electrode  surface  in  cm2. 

The  kinetic  electrode  parameters  aa,  aCt  Jo-a  and  Jo-o  have  a 
value-dependence  determined  by  the  material  choice  of  the  cata¬ 
lyst,  the  operating  temperature  and  the  electrolyte  type  and 
weight-weight  percentage  (wt.%).  In  the  HRI’s  alkaline  electrolyser 
(described  in  Table  1 ),  Nickel  (Ni)  at  99.99%  of  purity  has  been  used 
as  the  catalyst  for  the  cathode  and  anode  as  well  as  aqueous  po¬ 
tassium  hydroxide  (KOH)  at  30  wt.%  (optimal  weight-weight  per¬ 
centage  for  KOH  electrolyte  [7]).  The  approximations  of  Jo-c  and  Jo-a 
for  the  HRI  electrolyser  (data  reported  in  Refs.  [17—20])  can  be 
modeled  by  (7)  and  (8),  respectively. 

Jo-c  =  13.72491  -  0.090557  +  0.0905572  (7) 

a  =  30.4  -  0.206T  +  0.0003572  (8) 

The  expressions  describing  the  evolution  of  ac  and  aa  have  been 
found  using  experimental  data  published  in  Refs.  [17-20  for  the 


Table  1 

HRI’s  Electrolyser  parameters. 


Operating  information 

Description 

Max.  operating  voltage 

48-56  V 

Operating  pressure  (P) 

1  bar 

Operating  temperature  (T) 

0-353.15  K 

Hydrogen  production  rate 

1  Nm3  h-1  at  80  °C 

Electrical  power  reference 

5  kW  at  60  °C 

Electrolyser’s  component 

Material 

Membrane 

Zirfon 

Anode 

Nickel  99.99% 

Cathode 

Nickel  99.99% 

Electrolyte 

KOH  at  30  wt.  % 

Symbol 

Description 

Value 

Hcell 

Number  of  electrolytic  cells 

24 

dac 

Anode— cathode  gap 

2.5  mm 

dam 

Anode— membrane  gap 

1.25  mm 

dcm 

Cathode— membrane  gap 

1.25  mm 

Sm 

Membrane  surface  area 

0.03  m2 

Membrane  thickness 

0.5  mm 

Sa>  Sc 

Anode  and  cathode  surface  area 

0.03  m2 

ea,  ec 

Anode  and  cathode  thickness 

2  mm 

G,  lc 

Anode  and  cathode  height 

45  cm 

Jo-c 

Cathode  exchange  current  density 

(7) 

Jo-a 

Anode  exchange  current  density 

(8) 

ac 

Cathode  transfer  coefficient 

(9) 

aa 

Anode  transfer  coefficient 

(10) 

60 
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given  catalyst  materials  and  electrolyte;  the  correlations  used  for 
the  HRI  electrolyser  are  given  by  (9)  and  (10)  [7]. 

ac  =  0.1175  +  0.00095  T  (9) 

aa  =  0.0675  +  0.00095  T  (10) 

The  bubble  effect  on  the  electrode  surfaces  is  included  in  the 
electrode  activation  overvoltage  definition. 

The  activation  overvoltage  of  the  electrodes  taking  into  account 
6  is  given  by  (11). 

Vact-0  =  «<*(^)  (11) 

The  activation  overvoltage  of  the  electrodes  can  be  separated 
into  two  different  parts  as  shown  in  (12). 


ffK0H-e  =  n  -  e)\  (16) 

^KOH-free 

where  0KOH-e  is  the  new  ionic  conductivity  of  the  electrolyte  as  the 
result  of  the  bubbles  presence  in  the  bulk  electrolyte.  The  rela¬ 
tionship  between  e  and  6  under  specific  conditions  as  reported  in 
Ref.  [24]  is  given  by  (17). 

«  =  \0  (17) 

Re ie  is  defined  as  the  sum  of  the  bubble-free  electrolyte’s  resis¬ 
tance  (Reie-free)  and  the  resistance  due  to  bubbles  in  the  electrolyte 
(Re ie-,),  (18). 

Rele  —  ^ele-free  +  ^ele-e  (18) 

The  electrical  resistance  formula  is  also  applied  to  calculate  the 
resistance  due  to  the  electrolyte  (Reie-free)»  as  shown  in  (19). 


Vact-J 


^act  +  Vq 


+  Mog(l  -  0) 


(12) 


Me- free  — 


1 

°KOH— free 


-tam 

sT 


ucm 


(19) 


where  the  first  part  (Vact)  is  related  to  the  overvoltage  as  the  result 
of  the  activation  of  the  electrodes,  and  the  second  part  (V$)  repre¬ 
sents  the  bubble  overvoltage  due  to  the  coverage  of  electrode 
surface. 


where  dam  and  dcm  are  the  distances  between  the  electrodes  (anode 
and  cathode,  respectively)  and  the  membranes. 

Finally  the  definition  of  Re \e.£  is  given  by  (20). 


2.1.3.  Electrical  resistive  phenomenon 

Each  component  in  the  electrolytic  cells  is  modeled  as  an  elec¬ 
trical  resistance.  Then,  Ra  and  Rc  are  the  anode  and  the  cathode 
resistances;  Rkoh  is  the  resistance  of  the  electrolyte  (in  this  case 
KOH),  which  is  estimated  by  correlations  taking  into  account 
temperature  and  concentration;  Rmem  is  the  resistance  of  the 
membrane  separator.  A  bubble  phenomenon  affecting  KOH  re¬ 
sistivity  is  included  in  this  approach. 

2.13.1.  Electrodes.  The  electrical  conductivity  of  the  Ni  (oni),  in 
S  cm-1,  is  defined  by  (13)  as  a  function  of  T. 

aNi  =  6000000  -  2796507  +  53272  -  0.3805773  (13) 

The  resistances  of  the  electrodes,  Ra  and  Rc,  as  a  function  of  their 
geometry  and  conductivity  are  given  by  (14). 


Ra 


±(j*),Kcm'(j<) 

°Ni  V5a/  °Ni  V5c/ 


(14) 


where  Sa  and  Sc  are  respectively  the  anode  and  cathode  cross- 
sections.  La  and  Lc  are  the  electrode  thicknesses. 


2.13.2.  Electrolyte.  The  ionic  conductivity  of  the  KOH  (aKOH-free)  in 
S  cm-1  depends  on  the  temperature  and  the  molar  concentration 
(given  by  the  m  in  mol  m-3)  of  the  electrolyte.  Empirical  charac¬ 
terizations  have  been  reported  in  Refs.  [21,22].  The  correlation 
published  by  Gilliam  et  al.  [21]  has  been  used  in  this  work.  This 
correlation  is  shown  by  (15). 


^ele-e  —  ^ele-free  ^/2  ^  (20) 

The  existence  of  Re \e.£  depends  on  the  presence  of  a  current 
passing  through  the  alkaline  cells  and  is  calculated  using  Reie-free 
and  e  as  shown  in  (20).  If  the  operating  current  is  equal  to  zero,  the 
Reie-f?  is  equal  to  zero  too. 

2.1.4.  Membrane 

The  HRI  electrolyser  uses  a  non-organic  membrane  based  on 
Zirfon.  An  empirical  relation  between  electrical  resistance  and 
temperature  for  a  0.5  mm  thickness  Zirfon  membrane  has  been 
reported  in  Ref.  [25]  in  the  form  of  (21 ),  where  Sm  is  the  membrane 
surface  in  cm2. 

_  0.060  +  80er/50 

mem  -  lOOOOSm  (  j 


2.2.  Physics  model  validation 

The  physics  model  has  been  developed  based  on  a  “mean” 
model.  This  model  simulates  only  one  electrolytic  cell,  and  sup¬ 
poses  that  all  the  cells  contained  in  the  electrolyser  (nceii)  have  the 
same  physical  performance  and  behavior.  The  relationship  between 
1/ceii  and  the  electrolyser’s  operating  voltage  Veie  is  thus  given  by 
(22). 


Vele  —  ^cellKiell 


(22) 


<7 KOH-free  =  -2.04m  -  0.0028m2  +  0.005332m7  +  207.2  ™ 

+  0.001043m3  -  0.0000003m2T2 

(15) 

The  gas  bubbles  created  on  the  electrodes’  surfaces  go  into  the 
bulk  electrolyte  affecting  the  OKOH-free-  This  effect  can  be  modeled 
by  using  the  Bruggman  equation  (16)  reported  in  23]. 


In  this  physics  model,  all  the  information  needed  to  set  up  the 
electrolyser  cell  model  is  divided  into  two  parts:  the  first  part  is 
devoted  to  the  operating  parameters,  and  the  second  part  contains 
all  the  information  relative  to  the  electrode  material,  membrane, 
electrolyte  concentration  and  cell  geometry. 

An  atmospheric  electrolyser  has  been  used  to  make  a  first 
validation  of  the  physics  model.  Table  1  presents  a  general 
description  of  the  HRI’s  alkaline  electrolyser  used  in  this  validation. 
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Fig.  1.  Electrolyser’s  model  validation  at  different  temperatures  for  the  HRI’s  electro- 
lyser  at  atmospheric  pressure. 


Fig.  1  presents  the  comparison  between  the  experimental  and 
the  simulation  data  of  i—V  curves  at  different  temperatures  for  the 
HRI  electrolyser. 

Experimental  data  measured  at  HRI  have  been  used  to  validate 
the  proposed  physics  model.  The  i—Ve ie  simulation  results  shown  in 
Fig.  1  for  different  operational  temperatures  have  been  calculated 
with  the  equations  presented  so  far. 

For  purposes  of  comparison  to  our  earlier  work  [3],  Fig.  2  shows 
the  same  experimental  data  as  in  Fig.  1  but  fitted  with  the  pa¬ 
rameters  of  that  older  model. 

The  results  using  the  physics  model  are  in  better  agreement 
with  the  experimental  data  than  those  using  the  model  proposed  in 
Ref.  [3].  The  accuracy  of  the  new  results  is  better  irrespective  of  the 
operating  temperature.  The  physics  model  is  able  to  describe  the  i— 
Veie  response  of  the  HRI  electrolyser  from  the  start-up  period  (low 
current)  to  maximal  operating  current  at  different  temperatures, 
which  means  that  it  is  an  improvement  of  the  way  of  modeling 
alkaline  electrolysers  operating  at  atmospheric  pressure. 

The  aim  of  this  comparison  is  to  demonstrate  the  versatile 
character  of  the  model  proposed  in  this  paper  and  its  ability  to  be 


Table  2 

Phoebus’s  Electrolyser  parameters. 


Operating  information 

Description 

Max.  operating  voltage 

30-40  V 

Operating  pressure  (P) 

7  bar 

Operating  temperature  (T) 

0-353.15  K 

Electrical  power 

5  kW-26  kW 

Electrolyser’s  component  Material 

Membrane 

NiO  diaphragms 

Anode 

Ni/Co304/Fe  on  performed  Ni-plate 

Cathode 

C— Pt  on  performed  Ni-plate 

Electrolyte 

KOH  at  30  wt.% 

Symbol 

Description 

Value 

ftcell 

Number  of  electrolytic  cells 

21 

dac 

Anode— cathode  gap 

Zero  gap 

dam 

Anode— membrane  gap 

Zero  gap 

dcm 

Cathode— membrane  gap 

Zero  gap 

Sm 

Membrane  surface  area 

0.25  m2 

used  to  simulate  the  behavior  of  electrolysers  even  at  different 
operating  temperatures. 

The  electrolyser  model  presented  by  Ulleberg  in  Ref.  [2]  has 
been  used  widely  in  the  literature  to  simulate  electrolysers  in  any 
kind  of  simulation  systems.  That  model  is  based  on  a  semi- 
empirical  approach  and  is  however  restricted  to  a  specific  alka¬ 
line  electrolyser.  The  Phoebus  electrolyser  [2]  has  been  used  to 
validate  the  model  proposed  by  Ulleberg.  A  second  validation  for 
our  physics  model  has  thus  been  done  using  the  Phoebus  electro¬ 
lyser.  Table  2  presents  a  general  description  of  that  alkaline 
electrolyser. 

The  interest  of  simulating  the  Phoebus  electrolyser  using  the 
physics  model  proposed  in  this  paper  is  to  show  the  response  of  the 
physics  approach  using  a  pressurized  electrolyser  and  a  zero  gap 
technology  in  comparison  to  the  results  with  the  model  proposed 
in  Ref.  [2]. 

Simulations  for  the  Phoebus  electrolyser  have  been  carried  out 
using  physics  model.  Fig.  3  presents  the  simulation  data  with  the 
experimental  data  reported  in  Ref.  [2]  for  this  electrolyser. 

The  simulated  results  in  Fig.  3  show  a  good  correlation  with  the 
experimental  data.  Thus,  three  operating  temperatures  for  this 
electrolyser  have  been  chosen  in  order  to  compare  the  results  for 


Fig.  2.  Experimental  data  compared  with  the  model  propose  in  Ref.  [3]  for  the  HRI’s 
electrolyser  at  atmospheric  pressure. 


Fig.  3.  Electrolyser’s  model  validation  at  different  temperatures  for  the  Phoebus 
electrolyser  at  7  bar. 
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Fig.  4.  Comparison  between  Uileberg  approach  [2]  and  HRI’s  model  for  Phoebus 
electrolyser.  Vcen  for  different  temperatures  at  7  bar. 

the  Phoebus  electrolyser  using  physics  model,  Ulleberg’s  model  [2] 
and  experimental  data.  Fig.  4  presents  the  comparison  at  30  °C, 
50  °C  and  80  °C. 

As  shown  in  Fig.  4,  the  Uileberg  approach  and  the  physics  model 
are  equally  good  to  emulate  the  electrical  response  of  the  Phoebus 
electrolyser.  The  evolution  of  these  polarization  curves  for  all  three 
temperatures  are  in  a  good  agreement  with  the  experimental  data. 
The  difference  between  the  two  models  lies  in  the  fact  that  the 
physics  model  permits  tracking  the  evolution  of  the  pressure  P  and 
the  molar  concentration  m  with  changing  temperatures  while 
Ulleberg’s  model  considers  these  constant.  Since  in  real  electro- 
lysers  P  and  m  (and  maybe  other  “configuration  parameters”  as 
well)  change  with  T,  it  follows  that  our  physics  model  should  be  the 
better  model  to  predict  the  behavior  of  such  electrolytes. 

To  obtain  an  even  more  realistic  simulation  performance  from 
the  physics  model,  modules  involving  external  thermal  and  pres¬ 
sure  calculations  could  be  added.  In  this  way,  and  applied  to 
hydrogen  production,  the  efficiency  calculations  will  be  even  more 
exact. 

3.  Electrolyser’s  electrical  analogy  and  integration  with 
power  electronics 

Power  electronics  interfaces  are  used  to  regulate  the  power 
transfer  between  energy  sources  and  storage  media  to  the  elec¬ 
trolyser  [9,26].  An  electrolyser’s  electrical  analogy  has  been 
developed  in  order  to  create  an  electrically  programmed  interface 
between  the  power  electronics  systems  and  the  physics  model.  This 
interface  permits  to  couple  the  physics  model  to  the  electrical 
simulation  systems.  This  electrical  analogy  of  electrolyser  can  be 
useful  for  many  power  control  system  based  on  hydrogen 
production. 

3.1.  Electrical  analogy 

The  V—i  relationship  of  the  physical  effects  presented  and 
modeled  in  Section  2  of  this  work,  have  been  analyzed  in  order  to 
establish  similarities  with  the  electrical  behavior  of  some  electronic 
components.  One  simulation  has  been  carried  out  at  80  °C  and  1  bar 
using  the  HRI  electrolyser,  and  the  results  are  shown  in  Fig.  5.  These 
five  electronic  components,  programmed  using  MATLAB/Simulink/ 


SimPowerSystems®,  are  used  to  generate  a  circuit  representing  the 
electrical  performance  of  the  electrolyte  (see  Fig.  6). 

3.1.1.  Zener  diode 

The  V—i  characteristic  of  a  Zener  diode  in  reverse  bias  has  some 
similarities  to  the  voltage  required  by  the  electrochemical  reaction 
of  the  electrolysis.  The  Zener  diode  voltage  (BV)  equivalent  is 
described  by  (23). 

BV  =  -  (£th(T,  p,m)  +  V«_a  +  Ve_c)  (23) 

where  Eth  is  the  onset  breakdown  voltage;  V#_a  and  Vq-c  are  the  non¬ 
linear  terms  simulating  the  overvoltage  due  to  the  gas  bubbles 
adhering  to  the  gas-evolving  electrodes.  This  required  breakdown 
voltage  has  to  be  reached  before  current  starts  to  flow  through  the 
cells  to  initiate  the  electrolysis. 

3.1.2.  Rectifier  diodes 

The  activation  overvoltage  of  electrodes  is  similar  to  the  forward 
voltage  effect  of  a  general-purpose  (rectifier)  diode  in  forward  bias. 
The  forward  voltage  has  the  same  electrical  behavior  as  the  acti¬ 
vation  of  the  electrodes  in  an  alkaline  cell.  The  forward  voltage  is 
given  by  (24)  for  the  cathode-diode  and  by  (25)  for  the  anode-diode 
analogy. 
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^diode-c 

—  ^act-c 

(24) 

^diode-a 

=  Vact-a 

(25) 

3. 1.3.  Resistances 

Linear  resistive  physical  effects,  which  are  part  of  the  electro- 
lyser,  are  grouped  into  a  single  resistor  (. Rq )  as  shown  in  (26).  Rq  has 
a  linear  behavior  and  its  value  depends  only  on  T,  the  geometrical 
parameters,  and  the  materials  used  to  build  the  electrolytic  cells. 

Rq  —  +  ^mem  +  ^ele-free  (26) 

The  Re ie_e  is  used  as  the  non-linear  resistor  component  in  the 
electrolyser’s  electrical  analogy. 

3.2.  Integration:  DC/DC  converter  and  electrolyser’s  electrical 
analogy 


this  current  density  condition  is  reached  the  first  diode  is  added  to 
the  circuit.  Therefore,  current  is  defined  by  (28). 

i  =  -  BV  -  ydiode-c  (28) 

Rq 

The  activation  of  the  anode  is  determined  by  J  >  J0_ a.  Once  the 
current  density  satisfies  this  condition,  the  anode  diode  is  activated, 
and  hydrogen  and  oxygen  can  begin  to  be  produced.  The  current 
then  becomes  (29). 

j  _  VC  ~  BV  —  ^diode-c  ~  ^diode-a  (29) 

Rq 

Gas  production  within  the  cells  causes  bubbles  to  form  on  the 
electrodes.  Some  of  these  bubbles  remain  adhered  to  the  electrodes 
while  others  will/can  form  a  bubble  curtain  in  the  cross-section  of 
the  electrolyte.  This  gas  production  modifies  the  current  somewhat. 
Then  i  can  be  redefined  by  (30). 


The  electrolyser’s  cells  require  DC  for  operating;  in  this  work  a 
DC  bus  supplies  this  energy.  In  order  to  analyze  and  integrate  the 
electrolyser  model,  a  DC/DC  power  electronics  converter  has  been 
used  as  power  interface  between  DC  bus  and  electrolyser,  which 
allows  the  control  of  the  power  transfer  from  the  DC  source  to  the 
electrolyser  by  means  of  an  adjustable  duty  cycle  (a).  In  this  work,  a 
multiphase  buck  converter  developed  by  IRH  has  been  pro¬ 
grammed  using  MATLAB/Simulink/SimPo we rSy stems®  for  the 
analysis  of  the  transient  response  (see  Section  4.1),  and  an  ideal 
buck  converter  for  the  steady-state  analysis  (see  Section  4.2).  The 
electrical  analogy  of  the  electrolyser  connected  directly  to  the  buck 
converter  is  shown  in  Fig.  6. 

The  initial  state  of  the  electrolyser  is  defined  by  Rq  and  the  first 
term  of  BV  (simple  resistance  in  series  with  a  Zener  diode).  If  the 
value  of  Vc  (voltage  of  the  output  capacitor  of  the  buck  converter)  at 
the  start  is  lower  than  BV,  the  current  value  is  i  =  0.  Once  Vc  has 
exceeded  BV  a  current  i  appears,  therefore,  i  is  calculated  using  (27). 


VC-BV 


(27) 


The  i  is  measured  and  used  by  the  electrolyser’s  model  to 
calculate  all  the  new  values  (voltage  and  resistance,  as  shown  in 
Fig.  7)  of  the  electrical  components  of  the  electrolyser’s  electrical 
analogy. 

In  order  to  activate  the  electrodes,  the  current  density  J  must  be 
greater  than  the  values  Jo-a  and  Jo-c  respectively  for  each  electrode. 
The  activation  of  the  cathode  is  achieved  only  when  J  >  Jo-o  When 


I  _  Vc  BV  ^diode-c  ^diode-a  ^Q) 

Rq  +  ^ele-e 

The  Alkaline  Electrolyser  Simulation  Tool  ( AEST)  is  composed  by 
the  electrolyser’s  electrical  analogy  and  the  physics  model  as 
shown  in  Fig.  7.  In  AEST  the  i  value  is  recalculated  at  every  simu¬ 
lation  step  through  the  internal  feedback  calculation  loop  (/measured 
-  see  Fig.  7).  This  loop  measures  the  current  flowing  through  the 
AEST  and  controlled  by  the  DC/DC  converter;  this  is  used  to  esti¬ 
mate  the  new  values  (voltage  and  resistance)  of  the  electrolyser’s 
electrical  analogy  from  the  physics  model.  The  i  loop  guarantees  a 
good  calculation  of  the  electrolyser’s  operation  parameters  and  a 
correct  information  feed-back  in  the  coupling  interface. 

It  should  be  noted  that  this  electrolyser’s  current  loop  has 
been  used  in  Ref.  [8  to  control  an  electrolyser  in  a  wind-FI2 
hybrid  network.  The  model  used  in  Ref.  [8]  to  emulate  the  elec¬ 
trolyser  is  a  simplification  of  the  model  proposed  in  Ref.  [2].  This 
simplification  is  based  on  a  constant  voltage  source  and  constant 
resistance.  These  elements  were  used  to  perform  the  electrical 
analysis  of  the  electrolyser  current  calculation.  In  this  paper, 
the  AEST  is  a  more  realistic  module  for  the  electrolyser  emula¬ 
tion;  the  AEST  is  able  to  perform  in  electrical  models.  The  pro¬ 
posed  AEST  has  been  integrated  into  an  electrical  simulation 
system  using  an  electrical  analogy  with  a  closed  current  loop 
calculation  implemented  in  order  to  estimate  the  current  of  the 
electrolyser. 

The  next  section  of  this  paper  will  show  how  the  AEST  could  be 
used  to  analyze  the  start-up  phase  of  an  electrolyser  using  power 
electronics  interface  to  control  electrical  power,  this  kind  of  simu¬ 
lations  cannot  be  done  using  a  constant  model  of  the  electrolyser. 


Alkaline  Electrolyser  Simulation  Tool 
(AEST) 


Fig.  7.  The  Electrical  electrolyser’s  analogy  coupled  to  the  physics  model. 


4.  Simulations  results 

The  simulations  have  been  carried  out  using  the  AEST  and  the 
buck  converter  (as  shown  in  Fig.  6).  The  first  group  of  simulations 
have  been  done  in  order  to  evaluate  the  starting  up  phase  of  the 
electrolyser  and  the  influence  of  the  changing  of  V,  T,  and  P  on  the 
electrolyser’s  performance.  These  simulations  have  confirmed  the 
influence  of  the  operating  parameters  on  the  power  consumption 
and  on  the  operating  current  of  the  electrolyser.  The  start-up  period 
of  the  electrolyser  is  also  reported  in  this  section.  In  the  second 
group  of  simulations,  the  ideal  DC/DC  converter-AEST  has  been 
used  in  a  hypothetical  electrical  hybrid  (renewable  source  and 
diesel  generator)  topology  with  a  hydrogen  production  system  in 
order  to  observe  the  steady-state  response  of  the  system. 
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Table  3 

Multiphase  buck  converter’s  electrical  simulation  parameters. 


Symbol 

Description 

Value 

L\,  L.2,  E3  and  L4 

Inductance 

0.3  mH 

C 

Capacitor 

100  pF 

Vbus 

DC  bus  nominal  voltage 

48  V 

Freq 

IGBT  commutation  Frequency 

15  kHz 

N 

Number  of  cells 

4 

4.1.  Electrolyser’s  transient  response 

All  simulations  presented  here  have  been  made  using  the  HRI 
electrolyser  configuration  with  the  multiphase  buck  converter 
(described  in  Table  3)  used  in  Fig.  6.  Different  values  of  a  have  been 
applied  in  order  to  obtain  different  input  voltages.  The  current 
generated  by  the  imposition  of  a  voltage  ctl/bus  is  also  shown  in 
Fig.  8. 

The  impact  of  temperature  on  the  electrical  performance  of  the 
simulation  module  has  also  been  evaluated.  Several  simulations 
have  been  made  in  order  to  evaluate  the  electrical  response  of  the 
system  at  different  operating  temperatures  and  pressures.  The 
electrolyser  current— voltage  curve  for  one  cell  at  different  tem¬ 
peratures  is  shown  in  Fig.  9. 

These  results  confirm  that  the  temperature  influences  positively 
the  electrolyser’s  current  at  a  given  operating  voltage.  The  power 
absorbed  by  the  electrolyser  at  a  constant  temperature  of  80  °C  and 
different  operating  pressures  is  shown  in  Fig.  10.  When  the  oper¬ 
ating  pressure  is  increased,  the  current  drawn  by  the  electrolyser  is 
reduced  significantly. 


4.2.  Electrolyser’s  steady-state  response 

The  AEST  and  the  ideal  DC/DC  converter  have  been  integrated 
into  a  possible  hybrid  hydrogen  production  system  based  on 
renewable  energy  in  order  to  evaluate  its  steady-state  response. 
The  topology  of  the  theoretical  system  is  presented  in  Fig.  11. 

This  small  stand-alone  hybrid  system  includes  a  Wind-turbine 
generator  (WT)  and  a  Diesel  generator  (DG)  as  the  power  sour¬ 
ces;  where  the  DG  is  the  backbone  element  of  the  system.  This 
hybrid  system  represents  the  typical  case  in  stand-alone  areas, 


where  hydrogen  production  can  be  used  to  reduce  dependency  on 
fossil  fuels.  The  AC/DC  converter  is  supposed  to  be  an  ideal  con¬ 
verter.  The  electrolyser  is  integrated  into  the  system  as  a  power 
regulator  component  as  well  as  a  slave  dump  load.  Electrolysers 
have  recently  been  proposed  to  accomplish  this  challenge  in  small 
grids  [27-30],  using  the  surplus  energy  to  produce  hydrogen.  Load 
power  and  wind  speed  are  emulated  using  predefined  profiles. 

The  power  generated  by  the  HRI’s  wind  turbine  (Bergey  Excel 
10  kW)  is  determined  by  using  the  correlation  from  experimental 
data  (output  power  as  a  function  of  wind  speed),  as  presented  in 
Fig.  12. 

A  12  kW  DG  has  been  used  in  these  simulations;  in  order  to 
improve  the  DG’s  efficiency,  the  DG’s  operation  low  power  limit  has 
been  set  at  5  kW.  The  operating  conditions  for  the  DG  are  defined 
by  (31).  The  DG’s  power  demand  is  determined  by  (32). 

5  kW  <  Pdiesel  <  12  kW  (31) 

^diesel  =  A/vind  —  ^load  (32) 


Fig.  10.  Electrical  power  absorbed  by  the  electrolyser  at  different  operating  pressures 
and  at  80  °C. 
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Fig.  11.  Topology  of  hybrid  hydrogen  production  system  based  on  renewables. 


Steady-state  simulations  have  been  carried  out  in  MATLAB/ 
Simulink®  using  a  time  scale  of  1  s  (simulation)  equal  to  1  h  of  wind 
profile.  The  power  produced  by  the  WT  and  the  DG,  as  well  as  the 
power  demanded  by  the  load  for  400  h  of  system  operation  are 
shown  in  Fig.  13.  Under  these  conditions  the  DG  supplies  power 
constantly,  albeit  mostly  at  its  low  cutoff  value  of  5  kW,  and  higher 
only  occasionally  when  the  load  power  peaks  and  exceeds  what  the 
WT  can  produce  at  that  moment.  It  is  of  course  the  excess  power 
generated  by  the  WT  which  gets  channeled  to  the  hydrogen¬ 
generating  electrolyser  (or  elsewhere,  as  symbolized  by  the 
“dump  load”  in  Fig.  11). 

The  available  power  is  determined  by  (33)  and  shown  in  Fig.  14. 
^available  =  ^diesel  “F  f'wind  —  ^load  (33) 

This  amount  of  electrical  power  must  be  absorbed  by  a  regulator 
in  order  to  avoid  frequency  grid  perturbations  [12]. 

The  electrolyser  is  used  as  the  master  element  to  absorb  the 
power  excess  of  the  system.  However,  Pe ie  is  limited  by  T,  P  and  m. 
Hence,  a  dump  load  is  used  as  a  second  regulation  component  in 
slave  mode;  this  dump  load  should  absorb  the  amount  of  power, 
which  the  electrolyser  is  not  able  to  use  to  produce  hydrogen.  The 


electrolyser  and  dump  load  powers  must  be  equal  to  the  available 
power  as  shown  in  (34). 

^available  =  ^dump  +  Pe\e  (34) 

Peie  is  calculated  using  the  coupling-system  and  taking  into 
consideration  the  duty  cycle  a,  the  DC  bus  voltage  and  the  elec¬ 
trolyser  current  (Pe ie  =  a-Vbus‘0- 1  is  estimated  by  using  (30);  this 
means  that  the  electrolyser  power  in  its  steady-state  can  be  given 

by  (35). 


Peie  =  aV\ 


bus 


( fl^bus  BV  ^diode-c  ^diode-a 


RQ  +  Rele-£ 


(35) 


Peie-max  is  found  when  l/bus  —  Vf  and  P0ie-min  is  fixed  at  300  ^ J. 
The  electrolyser  operates  in  these  simulations  at  constant  T,  P  and 
m,  control  is  effected  using  different  values  of  a  as  per  the  following 
conditions. 


Case  1.  :  If  PaVaiiabie  <  Peie-min,  the  electrolyser  is  turned  off. 

Case  2.  !  If  Pele-min  <  Peie  <  Pele-max*  the  electrolyser  is  turned  on 

and  the  duty  cycle  used  to  control  the  buck  converter  is  computed 
using  (36). 


Fig.  12.  WT  characteristic  curve. 


Fig.  14.  ^available  to  be  absorbed. 


Power(W) 
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Fig.  15.  Peie  and  Pdump- 
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Fig.  16.  Vgie  at  different  temperatures. 


«  =  0.94^!™e  +  0.003 (Pavailable  -  Pele)  (36) 

1  ele 


Case  3.  !  If  ^available  >  fele-max>  the  electrolyser  is  turned  on  and  the 

duty  cycle  of  the  buck  converter  is  set  to  its  maximum.  In  this  case, 
the  dump  load  is  activated  and  its  power  is  computed  by  using  (34) 
in  order  to  keep  a  power  balance. 

The  electrolyser  temperature  has  been  set  to  50  °C.  The  power 
absorbed  by  the  electrolyser  and  the  dump  load  are  plotted  in 
Fig.  15.  Clearly  PaVaiiabie  is  absorbed  mostly  by  the  electrolyser,  and 
the  dump  load  comes  in  to  play  only  infrequently.  Peie-max  for 
operating  at  50  °C  is  equal  to  4627  W.  The  variations  of  the  elec¬ 
trolyser  voltage  and  current,  for  three  different  operating  temper¬ 
atures,  are  plotted  in  Figs.  16  and  17. 


5.  Conclusions 

A  model  for  alkaline  electrolysers  has  been  developed 
considering  a  physics-based  approach.  This  approach  takes  into 
account  all  quantifiable  physical  descriptors  of  the  various  phe¬ 
nomena  appearing  on  the  electrodes  or  the  membrane  or  in  the 
electrolyte-contact  resistances,  threshold  voltages,  bubble  for¬ 
mation,  etc.  The  operating  temperature,  pressure,  electrolyte 
concentration  and  current  have  been  defined  as  the  inputs  of  the 
electrolyser  model.  This  proposed  physics  model  is  able  to  repli¬ 
cate  the  i—V  characteristic  curve  as  well  as  the  electrolyser’s 
performance. 

An  electrical  analogy  of  the  electrolyser  consisting  of  series 
arrangement  of  resistances  and  diodes  has  also  been  developed, 
each  of  these  components  reflecting  the  various  physical  effects 
taking  place  inside  the  electrolyser.  The  electrolyser’s  electrical 
analogy  coupled  with  the  physics  model  have  been  integrated  into 
an  electrical  simulations  system  called  an  Alkaline  Electrolyser 
Simulation  Tool  (AEST).  The  integration  of  the  given  power  elec¬ 
tronics  devices  and  AEST  using  the  Flydrogen  Research  Institute 
electrolyser  configuration  is  explained  in  this  paper,  as  well  as  the 
i—V  characteristics  resulting  from  the  simulation.  Comparison  with 
experimental  values  is  excellent.  Using  this  new  approach  the 
current  i  becomes  a  reaction  of  the  process  (imposing  a  V  on  the 
electrolyser),  therefore  the  imposed  V  acts  as  an  input  in  the  AEST. 

Finally  the  proposed  AEST  has  been  used  to  emulate  the  physical 
dynamics  of  a  hybrid  electrical  system  consisting  of  a  diesel 
generator,  a  load  network,  and  the  electrolyser  used  to  produce 
hydrogen  during  those  phases  of  the  system  when  excess  power  is 
available.  The  AEST  is  able  to  predict  hydrogen  production  using  the 
simulated  operating  electrolyser’s  current,  to  emulate  the  impact  of 
the  electrolyser  in  power  electronics  and  the  impact  of  electrical 
sources,  power  electronics,  temperature,  pressure,  electrolyte 
concentration  on  the  physical  and  electrical  response  of  the 
electrolyser. 

AEST  is  a  useful  tool  for  predicting  and  emulating  the  alkaline 
electrolyser’s  performance  and  the  hydrogen  production  in  long¬ 
term  (steady-state)  scenarios  as  well.  It  could  be  used  as  an 
excellent  tool  to  evaluate  control  strategies  in  systems  with  wind 
energy  influences.  The  number  of  required  parameters  used  in  the 
model  permits  the  characterization  of  any  alkaline  electrolyser. 
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